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Cavity optomechanical systems are being widely developed for precision force and displacement measure¬ 
ments. For nanomechanical transducers, there is usually a trade-off between the frequency (/m) and quality 
factor (Qm), which limits temporal resolution and sensitivity. Here, we present a monolithic cavity optome¬ 
chanical transducer supporting both high /m and high Qm- By replacing the common doubly-clamped, 
813 X 4 nanobeam with a tuning fork geometry, we demonstrate devices with the fundamental /m ~ 29 MHz 
and Qm ~ 2.2x10^, corresponding to an fMQM product of 6.35x10^^ Hz, comparable to the highest values 
previously demonstrated for room temperature operation. This high /mQm product is partly achieved by 
engineering the stress of the tuning fork to be 3 times the residual film stress through clamp design, which 
results in an increase of /m up to 1.5 times. Simulations reveal that the tuning fork design simultaneously 
reduces the clamping, thermoelastic dissipation, and intrinsic material damping contributions to mechanical 
loss. This work may find application when both high temporal and force resolution are important, such as in 
compact sensors for atomic force microscopy. 


Cavity optomechanical systems are being developed 
for many applications in precision force and displace¬ 
ment measurements^’^. Monolithic systems in which 
nanomechanical transducers are combined with inte¬ 
grated optical readout have been developed in geome¬ 
tries where optical resonances and mechanical modes are 
CO- located within the same physical structure^”^, and in 
systems for which optical and mechanical modes are sup¬ 
ported by different physical structures and are near-field- 
coupled^’^. Such near-held coupling enables the mechan¬ 
ical resonator size to be scaled down to the nanoscale 
while maintaining high displacement sensitivity^’^in con¬ 
trast to far-held optical readout, where dihraction ehects 
limit the mechanical resonator size that can be sensitively 
detected^^. For a given desired mechanical stihness (de¬ 
termined by the force sensing application), a nanoscale 
cantilever can have much higher resonant frequency /m 
( and therefore transduction bandwidth/temporal resolu¬ 
tion) than a microscale counterpart, due to its smaller 
ehective motional mass (m). Such a high frequency me¬ 
chanical resonator would ideally exhibit a high mechan¬ 
ical quality factor (Qm), as the force sensitivity scales 

as 1 /(/m^Qm^) (Ref. 11). However, there is usually a 
tradeoh between /m and Qm because shrinking down 
the resonator size comes at the expense of a reduction 
in Qm, due to increased clamping losses^Here, we 
demonstrate an integrated silicon nitride cavity optome¬ 
chanical system where high /m and Qm are simultane¬ 
ously achieved. Using microdisk optical resonators with 
intrinsic optical quality factor Qo > 6 x 10 ^ for readout, 
we develop a doubly-clamped tuning fork geometry as 
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the mechanical resonator, where we take advantage of 
elastic wave interference to limit the mechanical loss, 
while retaining the high tensile stress characteristic of 
stoichiometric Si 3 N 4 . We investigate different clamp de¬ 
signs that can increase the tensile stress by up to 2.9 
times and hence tune /m on chip, as well as surface 
treatment to improve Qm- Devices with /m ~ 29 MHz 
and Qm ~ 2 . 2 x 10 ^ are presented, corresponding to an 
/mQm product of 6.35x10^^ Hz. 

Doubly-clamped silicon nitride nanobeams have been 
extensively explored as nano mechanical resonators, with 
the high residual tensile stress (~1 GPa) produced by 
low-pressure chemical vapor deposition (LPCVD) of sto¬ 
ichiometric silicon nitride (Si 3 N 4 ) on silicon enabling 
Qm > 10 ^ to be achieved for MHz frequency modes^^. 
In the context of cavity optomechanics, such high-QM 
structures have been evanescently coupled to silica mi¬ 
crodisk resonators^ and incorporated within silicon ni¬ 
tride nanophotonic circuits^^. Here, our goal is to de¬ 
velop a 813 X 4 nanobeam-microdisk optomechanical sys¬ 
tem in which the mechanical frequency /m is increased 
to the tens of MHz range while maintaining high Qm- 
To do so, we replace the commonly used single beam res¬ 
onator with a tuning fork mechanical resonator, as shown 
in Fig. 1(a). 

For nanomechanical resonators operating in vacuum, 
there are several sources of mechanical energy dissipa¬ 
tion: clamping losses (QcZamp), thermoelastic dissipation 
(TED) (Qted), and material losses (Qmat)- The to¬ 
tal mechanical quality factor Qm can then be written 
as 1 /Qm = ^IQclamp + ^/QtED + 1/Q mat T ^jQother 
(Ref. 11 ). Clamping losses occur when the elastic energy 
radiates into its support structures^TED is the en¬ 
ergy dissipation due to strain-induced heating and the 
resulting temperature gradients^^’^^, and is expected to 
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FIG. 1. (a) False-colored scanning electron micrograph 

(SEM) of a fabricated device. The orange region indicates 
the Si 3 N 4 microdisk optical resonator while the blue region 
indicates the mechanical resonator and supporting structures, 
(b) Working principle of the device. The colored tuning fork 
shows the calculated mode shape of the first-order, in-plane, 
out-of-phase mechanical mode. The colored disk shows the 
TEi ,43 whispering gallery optical mode, (c)-(f) SEM images 
of four different mechanical resonator designs: (c) single beam 
(S), (d) tuning fork (E), (e) tuning fork with neck (EN), (f) 
high-stress tuning fork (HS). (g) Schematic of the character¬ 
ization system. PC, polarization controller; PD, photodiode; 
DAQ, data acquisition; ESA, electrical spectrum analyzer. 


dominate Akhesier damping for devices with our geom¬ 
etry/aspect ratios (so that the latter is not considered 
further in our discussion) In addition, losses caused 
by localized defect states both on the surface and in the 
volume of the material cannot be neglected^^. Compared 
with the modes of single cantilever devices^’^’^^ the in¬ 
plane, out-of-phase mechanical mode of the tuning fork 
structure (mode shape shown in Fig. 1(b)) enables the 
total force and moment at the outer end of the clamps 
to be zero. Therefore, the motion of the beams is ef¬ 
fectively decoupled from the clamps, resulting in an ex¬ 
pected lower clamping losses. This effect, which has been 
used in a variety of mechanical structures ranging from 
clocks and musical instruments (the tuning fork was in¬ 
vented in 1711^^) to atomic force microscope sensors^^, 
can also be understood as destructive interference of elas¬ 
tic waves, and similar ideas have been implemented in 
double disk cavity optomechanical structures'^. In con¬ 
trast to those structures, the doubly clamped geometry 
we adopt (Fig. l(a)-(b)) maintains the high residual ten¬ 
sile stress of the Si 3 N 4 film, which is important for main¬ 
taining high /m and 

To verify these benefits of the design, we fabricated 
cavity optomechanical transducers with four types of me¬ 
chanical resonators: single beams (S), turning forks (F), 
and tuning forks with necks (FN), as shown in Fig. 1(c)- 
(e). The tuning fork with and without neck geometries 
are nearly identical, with only the specifics of the clamp¬ 
ing geometry slightly modified by inserting a neck region. 
Each mechanical resonator was coupled to a microdisk 


optical resonator that is 15 /im in diameter, and the sep¬ 
aration between the optical and mechanical resonators is 
150 nm. The devices were fabricated in a 250 nm stoichio¬ 
metric silicon nitride film with an internal tensile stress 
of 1.1 GPa, using electron-beam lithography and dry 
etching processes similar to our previous work on high-Q 
Si 3 N 4 microdisk optomechanical resonators^^. The me¬ 
chanical resonator beams are 150 nm wide, and the beam 
lengths are varied between 12 /im and 40 /im to investi¬ 
gate the dependence of device performance on the beam 
length. Stress tuning is an effective way to increase both 
/m and Qm, and stress tuning of chip-based devices was 
previously realized by substrate bending^^. In the tuning 
fork design, a stress level within the mechanical resonator 
above that of the deposited film was achieved by increas¬ 
ing the clamp length on one side of the fork, for example, 
by 100 /im as shown in Fig. 1(f). Due to the initial un¬ 
balanced tensile forces in the clamp and beam (details 
available in supplementary materials^^), the stress is re¬ 
distributed after undercut, such that it decreases in the 
wide suspended clamp and increases by up to 3 times^^ 
in the attached narrow beams of the tuning fork. To our 
knowledge, this is the first time stress tuning has been 
achieved by a design-enabled, on-chip approach. This 
stress tuning method provides more flexibility and better 
control of engineering /m of nano mechanical resonators. 

For each mechanical resonator geometry shown in 
Fig. l(c)-(f), devices with three different beam lengths 
(12 /im, 20 /im, and 40 jam) were fabricated and char¬ 
acterized in vacuum (0.13 Pa) in order to evaluate both 
their /m and Qm- The characterization setup is shown 
in Fig. 1(g), with more details available in the supple¬ 
mentary material^^. The typical optical quality factor 
Qo is > 10^, as shown in Fig. 2(a). The optomechanical 
coupling coefficient {goMl‘^'^) of our device is calculated 
to be 140 MHz/nm by perturbation theory^"^ using fi¬ 
nite element method (FEM) determined mode solutions 
(Fig. 1(b)), with details provided in the supplementary 
material^^. Figure 2(b) shows the measured mechanical 
spectra of devices with the same beam length (20 jam) 
but different geometries. We find that the tuning fork 



FIG. 2. (a) Typical optical spectrum of microdisk resonator. 
The intrinsic optical Q factor is 6.07x 10^ ± 6x 10^. The 95 % 
confidence interval range for Qo is determined by a nonlinear 
least squares fit to the data, (b) Mechanical spectra of a single 
beam, tuning fork, and high-stress tuning fork structures with 
20 /xm beam length. The 95 % confidence intervals for Qm 
from the curve fitting are typically ±3 %, and are summarized 
in Fig. 3(a)-(b). 
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with and without neck increase Qm by 12 times and /m 
by 1.4 times the values of a single beam structure, respec¬ 
tively. The high-stress tuning fork device with a longer 
clamp produces an even more pronounced increase in Qm 
and /m, by 13.5 times and 2.1 times the values of a single 
beam structure, respectively. The stress-tuning induced 
/m increase is 1.5 times compared with regular tuning 
forks. We note that this device exhibited the highest 
ImQm product (6.35x10^^ Hz) for any of the structures 
we have investigated. 

The thermomechanical force noise is also decreased us¬ 
ing the tuning fork geometry. Assuming the force noise 
is frequency independent, its spectral density can be es¬ 
timated by Sf = ^kkBT/{27TfMQM) (Ref. 11), where 
k is the cantilever effective spring constant, ks is the 
Boltzmann constant, and T is the temperature. For 
20 /im long single beam, fork, and high-stress fork de¬ 
vices, the FEM calculated k are 13.23 N/m, 26.27 N/m, 
and 60.58 N/m, respectively. Using the calculated k 
values and the measured /m and Qm values, the cor¬ 
responding S]/‘^ are 0.40 fN/Hz^/^, 0.14 fN/Hz^/^ and 
0.16 fN/Hz^/^, respectively. These results show that, 
although k is increased by the tuning fork design, the 
high fMQM product still reduces the force noise by a 
factor of ^2.9 compared with the single beam devices. 
We also note that the mechanical damping coefficient 
7 = k/(2'KmfMQM) and the corresponding thermody¬ 
namic force noise Sf = AkBTm^ do not further decrease 
by the high stress tuning fork, which introduces an addi¬ 
tional stress increase of 2.3 times compared with regular 
fork^^, as the increase in /mQm is offset by the increase 
in k. 

The full parametric study of single beam, tuning fork, 
and tuning fork with neck structures is shown in Fig. 3, 
where the labels PI, P2, and P3 indicate experimen¬ 
tal results from different batches of devices. As shown 
in Fig. 3(a), the experimentally measured /m values 
match the FEM simulation results, indicating the mea¬ 
sured modes were the in-plane, out-of-phase mechanical 
modes. Eigure 3(a) further shows that the measured /m 
is higher in the tuning fork device architectures compared 
to single-beam devices for all beam lengths. This con¬ 
firms that the tuning fork structure helps to decouple 
motion of the beams from the clamps, which effectively 
increases the clamp stiffness. 

According to Eig. 3(b), the tuning fork design helps to 
increase Qm for all the beam lengths in the experiment. 
We believe the reduction of clamping losses is one of the 
primary reasons for the experimentally measured Qm im¬ 
provements, based on the following discussion. Eirst, we 
plot measured Qm values for both the out-of-phase and 
in-phase mechanical modes for the tuning fork structure 
in Eig. 4(a). In all the cases, the out-of-phase modes 
have higher Qm values than the corresponding in-phase 
modes, which is expected because the in-phase modes do 
not produce destructive elastic wave interference in the 
clamping regions. Next, we note that the tuning fork de¬ 
sign is more effective in achieving increased Qm m devices 
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FIG. 3. (a) FEM calculated and experimentally measured 

frequencies of the mechanical resonators, where the labels 
PI, P2, and P3 indicate experimental results from different 
batches of devices. The horizontal axis labels S, F, and FN 
stand for single beam, tuning fork, and tuning fork with neck, 
respectively, while 12, 20, and 40 specify the beam length in 
micrometers, (b) Experimentally measured Q^^, FEM calcu¬ 
lated Q^ed^ FEM calculated Q^ed + Q^at- The 95 % 
confidence interval ranges for fM and Qm, determined by a 
nonlinear least squares fit to the data, are denoted by the 
solid black regions above each bar in the graphs, (c) FEM 
calculated normalized temperature gradient for tuning fork 
and single beam structures during vibration, (d) FEM calcu¬ 
lated elastic strain energy density on a logarithmic scale for 
tuning fork and single beam structures during vibration. 


with shorter beams. This is because the fraction of the 
modal mechanical energy within the clamping regions is 
comparatively smaller for longer beams, so that for long 
enough beams, clamping loss is likely no longer a domi¬ 
nant loss mechanism^^. We also fabricated devices with 
80 /im long beams (data not shown), with measured Qm 
comparable to the 40 /im devices. We note that very 
little difference is seen in the experiments between the 
/m and Qm values of the tuning forks with and without 
neck, indicating that the energy dissipation reduction is 
essentially the same for both. 

In addition to the clamping losses, we also consider 
whether the TED contribution to Qm may be improved 
by the tuning fork geometries, using EEM simulations 
in which material deformation is coupled to temperature 
gradients^^. Indeed, we do find that the TED contribu¬ 
tion in tuning fork devices is smaller than that in single¬ 
beam ones for all beam lengths, as shown in Eig. 3(b). 
Simulation results of the local temperature gradient are 
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FIG. 4. . (a) Experimentally measured and FEM calcu¬ 
lated out-of-phase (left) and in-phase (right) me¬ 

chanical modes. The FEM calculated mode shapes are also 
shown in the insets, (b) Experimentally measured before 
and after HE treatment. The 95 % conhdence interval ranges 
for /m and Qm , determined by a nonlinear least squares fit to 
the data, are denoted by the error bars in the graphs. Some 
of the error bars are smaller than the bar graph line widths. 


shown in Fig. 3(c). The localized deformation in the 
tuning fork structure helps to decrease the deformation- 
induced heating and cooling areas. Therefore, the heat 
flows driven by the temperature gradient become smaller 
in the tuning fork structure, which results in lower TED 
than that of the single-beam structure. In addition, 
calculation of TED for both the out-of-phase and in- 
phase modes in tuning fork geometry (Fig. 4(a)) confirms 
that the out-of-phase modes have lower TED. However, 
the TED improvements in simulation are much smaller 
than the total loss improvements observed in the exper¬ 
iment, which indicates that the Qm improvements are 
not mainly due to decreased TED. 

To simulate the effect of material losses, similar to 
a previously published work^^, we introduce an addi¬ 
tional imaginary part to the Si 3 N 4 Youngs modulus, so 
that it becomes E = Ei iE 2 ^ and determine the re¬ 
sulting quality factor Qmat from finite element simula¬ 
tions. Under the assumption that QIQdamp + ^/Qother) 
is small in the tuning fork geometries, we choose a value 
of E 2 = 20 MPa {El = 290 GPa) to achieve a rea¬ 
sonable agreement between experimental losses and sim¬ 
ulated {IIQted + 1/Qmat) for the tuning fork devices. 
This single choice of E 2 is consistent with the literature^^ 
and results in losses that do not exceed the experimen¬ 
tally measured total losses for any of the devices. The 
simulation results show that the tuning fork structure 
helps to decrease the damping due to the material defect 
losses {1/Qmat) i^ addition to decreasing the clamping 
losses. This can be explained by the simulated elastic 
energy density shown in Fig. 3(d), corresponding to the 
strain distribution during vibration. The deformation of 


the tuning fork structures is more localized in the center 
beam region, while the deformation in single beam struc¬ 
tures is distributed in both beams and clamps. With a 
smaller amount of deformed material, damping due to 
material loss is reduced in tuning fork structures. 

As shown in Fig. 3(b), differences between experimen¬ 
tally measured 1 /Qm and simulated (1 /Qt£;d+ 1 /Qmat) 
of single beam devices are much larger than those of the 
tuning fork devices, especially for shorter length beams, 
indicating higher clamping losses for single beam struc¬ 
tures. In total, we believe the clamping loss, damping due 
to material loss, and TED are all reduced in the tuning 
fork structure, through localized deformation resulting 
from elastic wave interference. 

From the characterization we also noticed that al¬ 
though they are nominally the same structures, sam¬ 
ples from chip P3 have lower Qm than those from other 
batches. We briefly consider whether this is due to sur¬ 
face losses, which are parts of the material losses. Such 
losses are caused by the friction processes resulting from 
surface defects and roughness^^’^^ and have been viewed 
as a potentially universal limiting loss mechanism in 
Si 3 N 4 nanomechanical resonators^^, might be reduced in 
our structures via surface treatment. We immersed the 
P3 devices into 1:10 hydrofluoric acid: water (HF:H 20 ) 
for 90 s to remove a thin layer of material that may have 
been damaged during the dry etching process. According 
to Fig. 4(b), Qm for all of the devices on this sample were 
improved by up to 2 times after the HF treatment, ap¬ 
proaching (but not exceeding) the highest values we have 
achieved in untreated samples (shown in Fig. 2). Further 
investigation is needed to determine whether a surface 
treatment can be applied to yield higher Qm values than 
the best untreated samples. 

In conclusion, we have developed silicon nitride tuning 
fork cavity optomechanical transducers. These structures 
simultaneously increase the resonant frequency (/m) and 
mechanical quality factor {Qm) of the fundamental me¬ 
chanical mode compared with single beam devices, by 
up to 1.4 times and 12 times, respectively, through an 
effective increase in the clamp stiffness and reduction in 
energy dissipation. By engineering the clamp geometry 
further, calculation indicates that we increase the ten¬ 
sile stress in the beams by 2.9 times compared with the 
intrinsic stress of the Si 3 N 4 film, resulting in an exper¬ 
imentally measured frequency increase of 1.5 compared 
with that of the regular tuning fork, while the mechani¬ 
cal damping remained unchanged. The highest measured 
ImQm product of 6.35x10^^ Hz at room temperature 
is on par with the highest values reported for doubly- 
clamped Si 3 N 4 beams. This tuning fork design with both 
high Qm and high /m can find applications in force sens¬ 
ing applications, including those that use active feedback 
cooling to damp mechanical motion^^’^^. 
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Supplemental Material 


I. FABRICATION 


The devices were fabricated (Figure SI) on commer¬ 
cially available Si 3 N 4 coated silicon wafers. The stoi¬ 
chiometric Si 3 N 4 layer was 250 nm thick and grown via 
low pressure chemical vapor deposition (LPCVD), with 
an internal tensile stress of ~ 1.1 GPa measured by the 
wafer bowing method. A layer of positive tone electron- 
beam (E-beam) resist was spin-coated on the Si 3 N 4 film, 
followed by E-beam lithography at 100 keV and develop¬ 
ment in hexyl acetate at 8°C. The patterns were trans¬ 
ferred to the Si 3 N 4 film by a CHE 3 /Ar /02 inductively- 
coupled plasma reactive ion etch (RIE). Note that this 
RIE step left a thin layer of Si02 on top of the exposed Si, 
preventing a successful KOH undercut. Thus, we carried 
out an additional SE 6 /C 4 E 8 inductively-coupled plasma 
RIE to remove the Si02 layer. The resist was then re¬ 
moved by Nanostrip. The fabrication process is finished 
by device undercut in a 20 % KOH bath at room tem¬ 
perature and a following N 2 blow-dry. 


(a) 


Stoichimetric LPCVD 



(b) E-beam lithography 

& CHF3/02/Ar dry etch 
- —E-beam resist 

_ ^ SiO^ 


(c) 


SF/C4H3 dry etch 



(d) 


KOH undercut 


:7x; 
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FIG. SI. Fabrication Process 
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II. CHARACTERIZATION SETUP 

Measurements of optical and mechanical modes were 
performed using the characterization system shown in 
Figure S2. Light from a 1425 nm to 1470 nm tunable laser 
was sent through a polarization controller and into an op¬ 
tical fiber taper waveguide^’^, which provided evanescent 
coupling into and out of the microdisk resonator. The 
device under test was housed in a vacuum chamber at 
a pressure ~ 0.13 Pa 10“^ torr). Motion of the me¬ 
chanical resonator introduced a time-varying effective re¬ 
fractive index of the microdisk and hence modulated the 
optical mode wavelength. Thus, fixing the laser wave¬ 
length onto the shoulder of the optical mode resulted in 
intensity-modulated light transmission from the optical 
resonator. This output light was split by a fiber coupler, 
where 10 % of the power was used for swept-wavelength 
spectroscopy of the optical cavity modes, while the re¬ 
maining 90 % was sent to a photodetector whose electri¬ 
cal output signal was spectrally analyzed. Each mechan¬ 
ical mode was manifested in a strong peak in the radio 
frequency (RE) spectrum. 


transmission 



mech, spectrum 


FIG. S2. Characterization system 

In this work, we used the fiber taper probe due to the 
flexibility it provided in enabling measurement of devices 
across multiple chips. The fiber taper was fabricated by 
heating and stretching an optical fiber to thin its diam¬ 
eter down to approximately 1 m. In addition, a local 
indentation (dimple) with 10 m radius of curvature was 
formed at the thinnest region of the taper so that it could 
be used for selective probing2. In future work focused on 
force sensing applications, it will be replaced by an inte¬ 
grated on-chip wave guide with pigtailed optical fibers. 

III. FINITE ELEMENT METHOD (FEM) SIMULATION 

In order to better understand the performance of the 
tuning fork design as well as the experimentally mea¬ 
sured spectrum data, we calculated the residual stress, 
mechanical modes, thermoelastic dissipation (TED), op¬ 
tical modes, and the resulting optomechanical coupling 
parameters in our devices in a commercial EEM software. 


The principles of the calculation, as well as the key re¬ 
sults of stress tuning, are described in detail below. 

A. Stress tuning by the clamp design 

All the devices are doubly clamped to maintain the in¬ 
trinsic high tensile stress in the LPCVD Si 3 N 4 film (1078 
MPa in our samples). The uniformly distributed stress 
in the film and in the device before undercut will be re¬ 
distributed after undercut. As a result, the free-standing 
clamps will shrink after the devices undercut, resulting 
in the thin beams being stretched, and hence the residual 
stress is increased (as shown in Eigure S3). The residual 
stress was calculated by a structural mechanics station¬ 
ary analysis, given a fixed boundary condition at the end 
of each clamp and the pre-stress in the Si3N4 film (1078 
MPa) as the initial condition. The result shows that for 
all the devices, the tensile stress in the beams became 
larger while that in the clamps became smaller. 

Eigure S3 shows the result of the tuning fork and sin¬ 
gle beam device with 20 fim beam length. Eor the single 
beam device, the stress is increased to 1316 MPa in the 
beam after undercut, while decreased to 197 MPa in the 
clamps. Similarly, the tuning fork has a stress of 1318 
MPa in the beam and 196 MPa in the clamps. The exag¬ 
gerated deformation shows the shrinkage of the clamps 
and the elongation of the beams. A further increase in 
the beam stress (and corresponding increase in mechan¬ 
ical frequency) can be achieved by modifying the clamp 
design. To that end, we also designed a high-stress fork 
which has a 100 /im long clamp on one side, as shown 
in Eig. 1(f) in the manuscript. In this device, the beam 
stress after undercut is calculated to be 3080 MPa, which 
is approximately 3 times the intrinsic Si 3 N 4 film stress. 
This stress increase is manifested by an experimentally 
measured /m increase of 1.5 times and Qm increase of 1.1 
times, respectively, compared with regular tuning forks 
(as shown in Eig. 2 in the manuscript; we note that the 
regular tuning forks themselves have the aforementioned 
1.3 times increase in stress compared to the deposited 
film). These results confirm that this on-chip, design- 
enabled stress tuning method can be used to widely ad¬ 
just the internal tensile stress and the resonant frequency 
of a nanomechanical resonator, which is of general impor¬ 
tance for most of nanomechanical resonators to achieve 
high /m and Qm- Einally, we note that these results are 
not specific to the tuning fork geometry, and that this 
stress tuning method is applicable to any doubly-clamped 
mechanical resonator with pre-distributed stress. 

B. Calculation of mechanical resonances 

The mechanical modes of the single beam and tuning 
fork resonators were calculated by solving the following 
eigenvalue problem for the equation of motion of the dis¬ 
placement held, with the residual tensile stress of the 
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0 1500 MPa 

FIG. S3. Calculated deformation and corresponding stress 
re-distribution for the tuning fork structure (a) before and 
(b) after undercut and those for the single beam (c) before 
and (d) after undercut. The initial beam length is 20 mum 
for all the devices shown. The color map shows the stress 
distribution. The deformation is exaggerated by 70 times. 

Si3N4 film and geometric nonlinearity taken into account: 
d^w 

^o^-V-<t = Fv (SI) 

Here, po is the initial material density, u is the dis¬ 
placement field, (T is the elastic part of second Piola- 
Kirchhoff stress tensor, and Fv is the volume force vector. 
A zero displacement boundary condition was applied to 
the outer ends of the fork clamps. 


The calculated Qted values are much larger than the 
experimentally-measured Qm, indicating that thermoe¬ 
lastic dissipation is not the dominant factor contributing 
to the observed energy decay rates. 

D. Calculation of optical resonances 

The whispering-gallery optical modes (WGMs) in the 
optical resonator can be described by the electromagnetic 
wave equation 

V X (V X E) = e{uj/c)^E (S3) 

Here, E denotes the electric field, e is the permittivity 
tensor, uj is the angular frequency, and c is the speed of 
light. 

Given a specific geometry, the electric field can be ob¬ 
tained by solving the eigenvalue problem in eq. (SI) 
around a specific wavelength. A perfectly-matched layer 
was used to simulate open boundaries. In addition, since 
most of the electric field is concentrated in the periph¬ 
ery of the resonator, the influence from the Si pedestal 
that supports the microdisk can be ignored in the calcu¬ 
lation. For a disk resonator with a diameter D ^ 15 fim 
and Si 3 N 4 thickness t ^ 250 nm, we found the TEi ^43 
WGM at 1439.89 nm, which has a good match with the 
experimental result (1441.79 nm). 


C. Calculation of TED of mechanical modes 


E. Calculation of the optomechanical coupling parameter 


The corresponding thermoelastic quality factor QTED 
for each mechanical mode was evaluated by the equation^ 

= v.(.vr) + <3-r,(|)^:| (S2) 

Here, Cp denotes the heat capacity of the solid at con¬ 
stant stress, T the absolute temperature, n tje thermal 
conductivity, Q the heat source per unit volume, Tq the 
reference temperature at which the strain and stresses 
take the initial values, and e the material stain tensor. 
Note the last term is a double dot product of two second- 
order tensors defined as T : tijUij, the result 

of which should be a scalar. The boundary condition 
we applied to the outer ends of the fork clamps is zero 
temperature deviation (isothermal boundary condition). 

The resonant frequencies of both the in-phase and 
out-of-phase in-plane modes show general agreement (to 
within 15 %) between simulation and experiment. The 
remaining discrepancy may be due to a number of fac¬ 
tors, including the fabricated beam widths not matching 
the design values used in simulation, and imprecise mod¬ 
eling of the clamp region and underlying silicon supports. 
The latter is particularly important for short single beam 
devices, for which the specifics of the clamp geometry 
has a significant influence of mechanical mode frequency. 


The optomechanical coupling parameter goM quanti¬ 
fies the shift in the optical resonance frequency uJq per 
unit displacement x of the mechanical resonator, that is, 
goM = duojdx. As presented in Ref.^, a perturbation 
theory for shifting material boundaries^ can be employed 
to derive the following expression for goM^ which can be 
calculated from the mechanical and optical modes as de¬ 
termined above: 


goM = 


f dA(Q ■ n) Ae|E|| p - A(e-i)|Dj 


2/dW|E( 


(S4) 


Here, Q is the normalized mechanical displacement 
field; n is the surface normal; E and D denote the 
modal electric and electric displacement fields, respec- 

tively, A = edielectric-<^air, and A(e-l) = dielectric- 

in which Cdieiectric ^ud Cair are the permittivities of the 
tuning fork and air, respectively. To calculate goM^ 
the optical and mechanical modes are separately calcu¬ 
lated, yielding E, D, and Q, and the integrals in eq. 
S4 are evaluated, with n obtained from the surfaces of 
the undeformed structure. The calculated /27r is 
140 MHz/nm for coupling between the first radial or¬ 
der, transverse electric (TE) polarized disk mode in the 





















1450 nm band and the fundamental, in-plane and out-of- 
phase mechanical mode of the 12 fim long tuning fork ge¬ 
ometry (/m ~ 33 MHz). The calculated goM values are 
only slightly smaller ~ 130 MHz/nm) for the 

longer 20 gm long tuning fork geometry (/m ~ 17 MHz). 

We also compare these result with the calculation of 
9om by a second approach, in which a deformable mesh 
is used to simulate the change in optical mode frequency 
with respect to mechanical displacement. Similar to 
our previous work^, we first simulated the relevant me¬ 
chanical mode of the tuning fork (in-plane and out-of- 
phase) to obtain its frequency and mode shape. Next, 
the tuning fork was deformed with the mechanical mode 
shape with a prescribed boundary displacement value. 
The optical modes in the disk with the deformed me¬ 
chanical resonator were then calculated by solving the 
electromagnetic eigenvalue problem. We then adjusted 
the boundary displacement value and re-solve the elec¬ 
tromagnetic eigenvalue problem. Repeating this proce¬ 


dure enabled qom to be obtained as the slope of the 
fitted frequency-modal displacement curve. The calcu¬ 
lated for the 12 /rm long tuning fork using this 

approach is 146 MHz/nm, which has a good match with 
that calculated by the perturbation theory approach. 
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